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Abstract. This work is concerned with a quantization of the Pais-Uhlenbeck oscillators from 
the point of view of their multi-Hamiltonian structures. It is shown that the 2n-th order 
oscillator with a simple spectrum is equivalent to the usual anisotropic n - dimensional oscillator. 
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1. This work is concerned with a quantization of the Pais-Uhlenbeck oscillators from 
the point of view of their multi-Hamiltonian structures. The family of such oscillators was 
introduced in the paper as a toy model to study field theories with higher derivative terms. 
Evolution of the Pais-Uhlenbeck 2n-th order oscillator is defined by the following equation 

i=l 

^ ■ where f2 = (co%), i = 1, ...,n is a set of positive parameters (frequencies). Equation (JTJ) can be 
obtained by variation of the Lagrangian 

(2) 

Introducing the natural oscillator coordinates 

*=n(^+". 2 >n(^+-?>* w 

i=l i=k+l 

and conjugate momenta Pais and Uhlenbeck proved that the Hamiltonian related to the 
Lagrangian (J2J) has the form 

1 n 
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It follows from this expression that the Hamiltonian is not positive defined. Recently, the 
quantization of the fourth order Pais-Uhlenbeck oscillator was curried out in some details by 
Mannheim and Davidson [2j (see [Sj, as well). The starting point for the authors of |2| was 
the Lagrangian (J2J) with n = 2. Using the Dirac method they obtained the Hamiltonian (see 
below first expression from (|12j) ) which is different in its form from the Pais-Uhlenbeck one 
(J2J. But the quantization of this new Hamiltonian leads to the Hilbert space H. containing 
negative norm states (this is because of improper sign of the commutator of a pair of creation 
and annihilation operators). This result is treated usually as the intrinsic property of theories 
with higher derivative terms. 

However, there are arguments that a satisfactory quantization of the Pais-Uhlenbeck 
oscillators can be curried out. For simplicity, let us consider the fourth order oscillator. In 
terms of the oscillator coordinates © 

d x o d x o 

and momenta (or velocities, which are identical to the momenta in this case) 
dqi dq 2 

Pl= dP P2= dF 

equation can be rewrited equivalently as a canonical system of Hamiltonian equations of 
motion for the two-dimensional anisotropic oscillator 



dqi 
dt 


_ dH c 
dpi 


= Pi, 


dpi ___ 
dt 


dH c __ 
dqi 


d<?2 
dt 


_ dH c 
dp2 


= P2, 


dp 2 __ 
dt 


dH c _ 
dp2 



-w?gi, 



(5) 



where 



He = ~(Pi + <"i<2?) + \{pI + "hi)- (6) 

The canonical quantization of this oscillator leads to the usual commutation relations among 
the creation and annihilation operators and to the standard Hilbert space with a positive norm. 
Thus, it is seems natural to deduce a relevant (for a quantization ) Hamiltonian formulation of 
the model directly from the equation of motion (jH) omitting a Lagrangian formulation. 

The difference of the Hamiltonians (jlj) (for n = 2) and (jBJ) tells us that the Pais-Uhlenbeck 
oscillator equation of motion can be obtained using nonequivalent Hamiltonian structures. In 
another words this oscillator is a bi-Hamiltonian system j3], In the case of the classical 
fourth order oscillator this fact was established in j^J. 

In the present paper we shall consider the quantization of the Pais-Uhlenbeck oscillators 
from the point of view of their multi-Hamiltonian nature. It will be shown that the 2n-th order 
oscillator with a simple spectrum (all frequencies from the set Q are different) is equivalent to 
the usual anisotropic n - dimensional oscillator. We shall start from the case of the fourth order 
oscillator and then the related formulas for the general case will be written. 
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2. The equation of motion of the forth order Pais-Uhlenbeck oscillator (n = 2 in (JTJ) has 
the form 

^ + K 2 + loD^ + ulu 2 2 x = 0. (7) 
It can be written in the form of a system of first order equations 

dXi d^2 dX3 dX4 . n 2\ 22 /\ 

where local coordinates of the "phase" space {x\ = x is the coordinate in 

the original space, x 2 is the velocity and so on). Integral curves of the vector field 

d 

V = x 2 di + x 3 d 2 + x A d 3 - ((w 2 + ul)x 3 + ulu^x 1 )d i , = — , (9) 

OXi 

are exactly solutions of the system (jHJ) (see, for instance, |Zj,|H])- The simplest way to obtain 
integrals of motion of the oscillator is to use the following equation 

V(H) = 0. (10) 

Note here that for this purpose in the paper |Hj was applied the known solution of Eq. (J2J). In 
the local coordinates Eq. (JTUj) has the form 

x 2 diH + x 3 d 2 H + x^diH — ((oJ[ + uj 2 )x 3 + uJlui\xi)d A H = 0. 

Since the components of the field V are homogeneous linear coordinate functions then the 
analytic solutions of Eq. (fTU|) are homogeneous polynomials in Xi. For the considered Pais- 
Uhlenbeck oscillator we have two independent positive defined quadratic integrals of motion 

#1 = Hx 4 + u\x 2 ) 2 + \uj\{x 3 + C^Xi) 2 , 

(11) 

H 2 = f (x 4 + ^fx 2 ) 2 + \u 2 (x 3 + ufxt) 2 . 

Taking into account the above definition of the oscillator coordinates qi and related momenta 
Pi one can see that the Hamiltonian (JHJ) is the sum of integrals Hi and H 2 from (jlljl 

He — Hi + H 2 , 

whereas the Pais-Uhlenbeck Hamiltonian is their difference 

Hpu = Hi — H 2 . 

Let us point out that the coordinates qi, q 2 and the integrals ((H)) are degenerate in the case 
uji = uj 2 . The simplest linear combinations of Hi and H 2 

r _ u\Hx-u)\H2 _ 1 2. .2~2 1 If. .2 1 , ,,2-S„2 _L 1-2 , ,2, ,2™ 

1 2 " (12) 

^ _ H 1 -H 2 _ 1, ,2, ,2„2 I If, ,2 _|_ , ,2\„2 1„2 _i_ „ „ 
^2 — ~ J2_J2 — 2 I 2 I + 2V 1 ~"~ UJ 2J X 2 — 2 3 ~"~ X 1 X 4) 
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gives us the pair of integrals which are distinct at u>i = u 2 , but C\ and C 2 are not positive 
defined. The integral C\ was obtained in [2] and it was taken as a Hamiltonian of the fourth 
order Pais-Uhlenbeck oscillator. The notation of [2] and ours are related by the following 
formulas 



X\ = p q , x 2 = julu>lq, x 3 = yuiu 2 x, x 4 

Now introduce a Poisson structure for the considered oscillator. Recall that a Poisson 
structure on a manifold is defined by a rank two contravariant tensor field II which is skew- 
symmetric 

U ij = _ u ji 
and satisfies the condition 

[[n, n]p : = ^(n™<9 m rp fc + u m3 d m u kl + ir k d m u ij ) = o. 

m 

Any Poisson structure induces a Poisson brackets {F, G} of arbitrary differentiable 
functions F, G on a manifold M. In local coordinates Xi the brackets are defined by the formula 

{F, G} = I&diFdjG. 

The vector field V is called locally Hamiltonian if there is a Poisson structure II such that 

£v(n) = (13) 

where £v defines the Lie derivative along the field V. In local coordinates the above relation 
has the form 

.<9IP 



dV l udV^ 



dx k dx k dx k 
where V i and II 4 - 7 are components of the vector field V and the Poisson tensor EL If there is 
such a differentiable function H that V can be represented in the form 

V H (.) = {•>#}, 

then H is called a Hamiltonian. In this case equations of motion take the canonical form 



§ = {*,*>• 



(14) 



Considering the relation (|13jl as an equation one can obtain a Poisson tensor n related to 
the field V. The simplest solution of this equation is a two-parameter nondegenerate Poisson 
tensor with constant components. Its components are represented in the following table (/ and 
g are arbitrary parameters) 






f 





9 


-/ 





-9 








9 





-u\ulf - (ujj + ul)g 


-9 





ufulf + (ul + ul)g 






(15) 
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This Poisson tensor Uf g is obviously skew-symmetric and satisfy the condition [[n,n]] l:, ' fc = 
in view of its constancy. It induces the following Poisson brackets for the coordinate functions 

{xi,x 2 }f, g = f, {X!, x 4 } /)S = g, 

(16) 

{x 2 , x 3 } fy9 = -g, {x 3 , x A } Lg = -u\u\f - {u)\ + uQg. 
It is not difficult to check that the dynamical equations (jHJ are generated by these brackets 

= {xi,H} fjg 

together with the Hamiltonian function 

H = aiH! + a 2 H 2 , (17) 

where the coeffitients have the form 

1 1 

ai ~ M^M4fT7y a2 - ~ (<4 - "DMf + 9)' (18) 

and can be choosing positive. Thus the dynamical equations (jHJ) (and the field V itself) are 
Hamiltonian ones and the two-parameter function H plays the role of a Hamiltonian. Remark 
that the integrals of motion C\ and C 2 are in involution in respect with these brackets 

{C 1 ,C 2 }to = 0. 

In the classical case the parameters / and g can be taken either arbitrary or fixed in any 
appropriate manner. For instance, we can put 

This choice gives the following nonzero components of the Poisson tensor Hi 

{x 1 ,X 2 }i = j-j, {x 3 ,X4}l = l. (19) 

UJ \U 2 

Other simple choice 

/ = 0, g = l, 

gives 

{xi,x 4 } 2 = l, {x 2 ,x 3 } 2 = -1, {x 3 ,Xi} 2 = -(uf + ul) (20) 
for the components of U 2 . Both mentioned Poisson structures generate the dynamical equations 

^i = { Xl ,C 1 } l = {x t ,C 2 } 2 . (21) 
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Thus the fourth order Pais- Uhlenbeck oscillator is a bi-hamiltonian system ( jl] , jH] , jH] ) with the 
Hamiltonians C±, C% and the Poisson structures III, n 2 . First from these structures has been 
used in j2j. Let us note in conclusion, that there is no a constant Poisson structure which 
generates the equations (JHJ) together with any of the Hamiltonians ([lljh 

3. The quasiclassic quantization of the Poisson structure (fi3j) can be considered as 
exact because of the dynamical equations (JHJ) of the fourth order Pais-Uhlenbeck oscillator 
are reducible to the canonically quantized form (jSJ) by the linear transformation. Assume that 
the hermitian operators Xi, i = 1, 4, related to the dynamical variables of the classical system 
Xi, subject the following commutation relations 



[xi,x 2 }f, g = ifrf, [xi,x 4 ]/, s = ihg, 

[x2, x 3 ]f, 9 = -ihg, [x 3 , x 4 ] ftg = -\hJ{ujlf - ih(uf + ufjg. 
Using these relations it is not difficult to show that the quantum dynamical equation 



(22) 



d£l . dx 2 dx 3 dXi 2 2 2 2~ 

~~dT = X2 ' ~~dT = XSl ~~dT = X4 ' -^ = -(^1+^2)^3-^1^2^! ( 23 ) 
can be represented in the Heisenberg form 

ih^=[x u H] f , g (24) 

where the quantum Hamiltonian H is obtained from the classical one (fTTjl by the replacement 
of the classical dynamical variables by the quantum dynamical variables X{ Xi, i = 1, ...,4. 
Remark, that there is no problems with the ordering of the quantum variables because all the 
terms of the form X{Xj in H include only commutative operators. It is easy to check that 
the quantum analogs of all the above considered classical integrals of motion Hi, i? 2 , C'i, C2 
commutate with each other. As in the classical case, fixing the parameters / and g one can 
obtain independent realizations of the quantum dynamical equations in the Heisenberg form 
with different Hamiltonians. For example, putting / = —^r^z, g = (see (fT5|) ) or / = 0, g = 1, 
(see f)2L)p). we obtain the equations 

[ ^~d~t = ^' ih ~d~t = ^' ^ 2 

respectively. In these equations for the calculation of the commutator [.,.]% it is necessary use 
the first fixed pair /, g, and for the calculation of the commutator [., .] 2 it is necessary use the 
second one. In both cases the roles of Hamiltonians play the operators 

Ci = -\uj\uj\x\ + + ul)xl + \x\ + oj\uj\xiX 3 , 

C 2 = \u\uj\x\ + \ + w^xl - \x\ + £2X4. 

Thus, we obtain that the quantum version of the fourth order Pais-Uhlenbeck oscillator is the 
bi-Hamiltonian system as well. 



Remarks on quantization of Pais- Uhlenbeck oscillators 7 

In view of the linearity of the quantum dynamical equations (|23j) one can easily writes out 
their operator solution 

x x = e~' 1Wlt a l + a 2 e~ iuJ2t a 2 + h.c. , 

x 2 = — iuJie~ lu)lt ai — iu 2 e~ luJ2i a 2 + h.c. , 

(26) 

x 3 = —uj\e ltJlt ai — u 2 e 1U2t a 2 + h.c. , 

x A = iujfe~ iullt a 1 + io;| e~ iu)2t a 2 + h.c. 

where we took into account the self-conjugacy of the dynamical variables Xj. Using the 
commutation relations ()22)1 we obtain nonzero commutators among the operators Oj, af, % — 1, 2 

h(ujf + g) , h(uff + g) 

[auail ~ ^^-'2^1-^1)- (27) 

^From the condition 

[ai,a+] = [a 2 ,a+] = 1, (28) 

imposed usually on commutators of creation and annihilation operators, we uniquely fix the 
parameters / and g 

f = l( UJl + u 2 ), g = -1^1 + ^1). (29) 

Substituting the solution (|27)|) in the expressions for Hi, H 2 , H, C\, C 2 and taking into account 
the commutation relations with defined above parameters ()29|) we obtain 

Hi = Ci + u 2 C 2 = 2u 2 (u 2 - ulf{atai + \), 

J ff 2 = C' 1 +^C' 2 = 2^(a; 2 2 -a; 1 2 ) 2 (a 2 + a 2 + i), 1 J 



and 



# = a 1 Hi + a 2 H 2 = hu 1 (afa 1 + ^) + ^2(^2" a 2 + -) (31) 
(7a = 2(a; 2 - u 2 ) (-ujf(afai + |) + w 2 4 (a+a 2 + ±) , 

(7 2 = 2(a; 2 2 -^)(^(a+a 1 + i)-^(a+a 2 + i)). 1 j 

These formulas show that in the case u>i 7^ u> 2 the quantum forth order Pais-Uhlenbeck oscillator 
is equivalent to the usual anisotropic harmonic oscillator. Hence, using the operators af, ai 
one can costruct the standard Hilbert state space 7i with the positive normalized basis vectors 

|Vw >= -yl=(air(atr\0 >, 
V m\n\ 

and the vacuum vector satisfying the condition 
adO >= a 2 |0 >= 0. 
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All quantum integrals of motion Q, Hi are diagonal in this basis. 

4. Let us return to the general case of the Pais-Uhlenbeck 2n-th order oscillator. Rewrite 
the equation in the form 

i=\ j=0 

where a™ is the j-th degree elementary symmetric polinomials in n variables uf,i=l,...,n 

a o=l, 0<j<n. 



l<ii<...<i,<n 



In these notation the equation is equalent to the following system of 2n first order differential 
equations 

^ = x i+1 , i = l,...,2n- 1, = -^o^^n-jj+i, xi = x. (33) 

3=1 

The system f!33p has n integrals of motion. In the Pais-Uhlenbeck variables qi, pi 

n—l n—1 

Qi = ^2 0-j~ 1 (i)x 2 ( n -j)-l, Pi = a j~ 1 ( i ) X *(n-j)> * = 1, (34) 
j=0 j=0 

where <7™ _1 (i) is the j-th degree elementary symmetric polinomials in n — 1 variables ul, k = 
1, i, ...,n (the variable oof is omited), these integrals of motion take the form of the harmonic 
oscillator energy 

In degenerate cases when some of the frequencies from the set Q coincide, related integrals Hi 
coincide too. Hence, as in the case of the forth order oscillator, the role of a Hamiltonian which 
generates the dynamical equations (J33|) must play an appropriate linear combination of H t . We 
put 

n 1 

H = / hHi, hi = — ™ — t— j 5V, (35) 

^n J -=u(4-4) 

where the factor with j = i in the product is omited. This Hamiltonian together with the 
Poisson structure n defined by the following nonzero components 

rr- l+1+2j = (-i)^-i +2 „ n = 2 £?=i w * 

(36) 

i = l,...,2n-l, 0< j< 

generates the dynamical equations (J33j) . In the formula (|36|) [a] denote an integral part of a 
number a. 
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The quantization of the Pais-Uhlenbeck 2n-th order oscillator we will realize by the above 
scheme. Assume that commutation relations among the operators Xi, i — 1, ...,2n, related to 
dynamical variables have the quasiclassical form 

[x h Xj] = i£IF. (37) 

Using these relations and the quantum Hamiltonian H obtained from the classical one (|33j) by 
the substitution x\ — > Xi, one can check that the quantum version of the dynamical equations 

— = X J+1 , l=l,...,2n-l J — = - 2^ ^ »2(n-i)+l 

j=l 

has the Heisenberg form (|24|) . Using the operator solution of these equations 

xi = ^e" iLJl 'ai + h.c, ^ = ^=i, i = l,...,2n, (38) 

and the relations (J37j) we obtain the commutation relations among the creation and annihilation 
operators a*, a/" 

[ai, at] = [at, dj] = [af, a+] = 0, i, j = 1, In. (39) 

The substitution of the solution of ()38|) in the Hamiltonian H gives us the Hamiltonian of 
n- dimensional anisotropic oscillator 

n 1 

# = /i^ i (a+a i + -). (40) 

i=l 

5. Let us conclude these remarks by some notes on the degenerate case of the fourth order 
oscillator. As it was pointed out above, under the condition uj\ = u 2 = uj it is convenient to 
exploit the independent integrals f!12D which take take form 



C s \ = \x\ + UJ 2 x\ + UJ^XiXz — ^uj 4 xI, 

C s2 = — \x\ + uj 2 x\ + ^uj^xj + x 2 x 4 . 



\x\ + u 2 x 2 2 + \u A x\ 1 

Substituting tux = uj 2 = in the formulas (J22J) we obtain the following commutation relations 
[xi,x 2 ] = i^/, [£1,^4] = i%, 

(42) 

[x 2 , x 3 ] = -1%, [x 3 , x 4 ] = -ih(u 4 f + 2u 2 g). 

Using these relations and the solution of the quantum dynamical equations ()23)1 in the 

degenerates case 

x t = e~' lwt a x + te~' lujt a 2 + h.c. , 

x 2 = -\ue~' luJt ai + (1 - iujt)e~ iuJt a 2 + h.c. , 

(43) 

x 3 = -w 2 e~ lwi *ai - u(2i + ujt)e' lujt a 2 + h.c. , 
x 4 = \u z e-' lUlt ai + cj 2 (-3 + itu)e~ iuJt a 2 + h.c. 
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we obtain the commutation relations among the creation and annihilaton operators eij, af 

+ h3u 2 f + g + + + ihu 2 f + g 

[ ai ,aJ\ = , [a 2 ,aj\ = 0, ai,aj = K,a 2 = = — . (44) 

4 CaT 4 u z 

Remark that the pair a 2 ,a\ commutes for any parameters /, g. It is useful to fix these 
parameters by the conditions 

[ai, af) — l, [a 2 , G&2 ] = [oi, aj] = [ a 2i a i1 = 0. 

We obtain g = —uj 2 f and f = The integrals of motions (|4*T|) in the terms of a^ 1- , have the 
form 

C s \ = 16w 4 a^"a2 + 4ia; 5 (a2a^ — aiaj ), 

C S 2 = — 8w 4 a^a2 — 4ia; 3 (a 2 a] f " — aiaj )■ 

Commutativity of the operators a 2 ,a^" tells us that constructing the state space of this Pais- 
Uhlenbeck oscillator it is necessary to take into account their classical character. 
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